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The exchange interaction between electrons and magnetic spins is consider- 
ably enhanced near interfaces, in magnetic multilayers. As a result, a dc current 
can be used to generate spin oscillations. We review theory and experimental 
evidence. The s-d exchange interaction causes a rapid precession of itinerant 
conduction-electron spins s around the localized spins S of magnetic electrons. 
This s precession has been observed directly (Weber et al., 2001) with electron 
beams through Fe, Co, and Ni films. Because of it, the time-averaged interac- 
tion torque between s and S vanishes. Thus, electrons do not interact at all 
with long-wavelength spin waves, in the bulk. An interface between a magnetic 
layer and a spacer causes a local coherence between the precession phases of 
different electrons, in a region within 10 nm from the interface (Slonczewski, 
1996; Berger, 1996). Also, a second magnetic layer with pinned S is used to 
"prepare" s in a specific direction, before the electrons cross over to the ac- 
tive magnetic layer. The current-induced drive torque of s on S in the active 
layer may be calculated from the spin current (Slonczewski) or from the spin 
imbalance A/i (Berger). This torque is equivalent to negative Gilbert damping, 
leading to instability of S. Spin current and A/i are proportional to each other 
(Berger, 2001) and can arise fom Fermi-surface translation, as well as from ex- 
pansion/contraction. In fields H normal to layers, the critical current I c for S 
instability is predicted to be proportional to the ferromagnetic-resonance fre- 
quency lu (consistent with the Tsoi et al. 1998, 2000, experiments). However, 
for in-plane H, due to elliptic S precession, I c is not proportional to to, but 
linear in H (Katine et al. 2000, experiments). Apart from the current-induced 
drive torque, an extra Gilbert damping is predicted near the interface even at 
zero current (Berger, 1996). It has been observed by ferromagnetic resonance 
(Urban et al., to be published). 

I. INTRODUCTION AND INTERACTIONS IN THE BULK 

Recently, it has been suggested that the exchange interaction between elec- 
trons and magnetic spins is considerably enhanced near interfaces, in magnetic 
multilayers. As a result, a dc current can be used to generate spin oscillations, or 
even to switch the spin direction from up to down. Using classical spin models 
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and simple mathematics, we start by describing the situation in bulk samples 
and single thin films. The effect of interfaces will be discussed in Sections II 
and III. Then, some of the existing experimental evidence will be covered. The 
analogy with a semiconductor laser is explained in Section VIII. 

In metallic ferromagnets, the s-d exchange interaction —2J s dS- S(f) couples 
1 the localized spin S (r) of 3d magnetic electrons to the itinerant spins s of 
4s conduction electrons. This interaction can be approximately represented by 
a classical field H s d = — J s d < S(r) > / acting on s. Here, [Ib is a Bohr 
magneton, and we have J s d — 0.1 eV and noH s d ~ 1000 T. 

We assume that a spin wave is present, i.e., that S is precessing around the 
field H (Fig. la), which includes the external field as well as demagnetizing and 
anisotropy fields. In the 1950s and 1960s, theoretical investigations of electron 
scattering by spin waves were often based on the first Born approximation. In 
that approximation, s-d exchange is assumed to be a small perturbation. In 
classical language, this means that s remains nearly fixed and nearly parallel or 
antiparallel to H (Fig. la), as S precesses around H . 

However, since H s d 3> H, s-d exchange is clearly too large for the Born 
approximation to be valid. A more correct picture 2 (Fig. lb) has s precessing 
rapidly around the instantaneous S direction, as S itself precesses more slowly 
around H. Depending on initial conditions, in many cases s follows closely 
the instantaneous S direction, namely, s is nearly parallel to S for a spin-up 
electron (Fig. lb) and antiparallel for spin-down. This could be called the 
strong-coupling (or adiabatic) approximation. Because the angle between s and 
S is now smaller (Fig. lb), the instantaneous value of the exchange torque 
between them is much reduced. Also, because of the rapid s precession around 
S, that torque changes sign frequently. Finally, different conduction electrons 
precess out of phase with each other in a bulk sample, so that their torques 
tend to cancel. We conclude 2 that electrons do not interact practically at all 
with spin waves through ordinary s-d exchange, at least for long-wavelength 
spin waves! 

This interesting conclusion must be modified 2 if the electron collision rate 
with solutes or phonons is high and comparable to the s precession frequency 
u> s . Also, if the spin- wave wavelength is short and comparable to vp/uj Sl where 
vp is the electron Fermi velocity; indeed, a current-induced exchange torque 
leading to a spin-wave instability has been derived 3 quantum-mechanically for 
finite spin-wave wavelength. The present paper will emphasize the limit of long 
wavelength (i.e., uniform S precession) and long electron collision time. 

The precession of s around S has been observed directly 4 by passing a 
polarized electron beam through a film of Fe, Co or Ni of thickness 0-10 nm 
(Fig. lc). The magnetic film is supported by a gold foil of thickness 20 nm. 
The electron beam is generated by a gun having a GaAs photocathode. The 
electron spin s is originally perpendicular to the fixed magnetic spins S of the 
uniformly magnetized film. After passing through the film, s is measured, and 
is found to be rotated by an angle e around S. The value of e is proportional to 
magnetic-film thickness. In cobalt, for example, the proportionality constant is 
19°/nm. In these experiments, the incident electron energy exceeds the Fermi 
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level in the film by only ~ 4 — 7 eV, and its exact value is found to have little 
effect on e. 

II. EFFECT OF REFERENCE LAYER AND INTERFACE 

S-d exchange torques exerted by conduction electrons on precessing magnetic 
spins S*2 in layer i*2 (Fig. 2a) can be increased 5 ' 6 over their bulk value, by 
introducing another magnetic layer F\ separated from F2 by a nonmagnetic 
spacer N. The magnetic spins S\ in F\ are assumed to have a fixed direction 
parallel to the precession axis of S2. If the F\ thickness is at least a local 
spin-diffusion length, this layer is thick enough to control the ratio of spin-up 
current to spin-down current in the multilayer, and to give to this ratio a value 
as different from one as possible (as needed 6 for current-driven experiments). 
As an example, the spin-diffusion length is 44 nm in Co nanowires at 77 K. 
A second role of layer F\ is to relax the spin s of electrons in F\ to a direction 
parallel (or antiparallel) to Si, under the combined influence of s-d exchange, 
spin-orbit interaction and random collisions. Thus, Fi prepares s in a fixed 
direction, and could be called the "reference" layer. 

Some electrons originating in Fi cross over to F2 (Fig. 2a). By solving the 
Schroedinger equation 5 > 6 > 8 , while assuming static and uniform S2, one finds that 
the s direction changes somewhat as the electron crosses the interfaces between 
these layers, but remains roughly parallel or antiparallel to S±. Thus, the angle 
between this s and S2 in F2 is large, and comparable to the 5*2 precession-cone 
angle 9 (Fig. 2a). In turn, this implies a large instantaneous torque between s 
and S2, as in the case of the Born approximation of Fig. la. But, as the electron 
propagates inside F2, away from the N — F2 interface, s precesses rapidly around 
§2 (Fig. 2b). As a result, the total torque exerted on all magnetic spins in F2, 
and obtained by summation over all s locations in F2 , is not in the tangential 
direction shown as f n in Fig. 2a, as one might naively think. Rather, it is 5,6 
in a direction Td radial to the 5*2 precession orbit; one exception being the case 
of a very thin layer F 2 discussed in Section VI. 

Electrons with opposite s produce opposite torques. Hence, after summation 
over all electrons, the total torque fd is found to be zero in equilibrium. An 
electric current I normal to layers is needed to populate certain electron states 
near the Fermi level and to depopulate others, resulting in a torque proportional 
to I. It is also approximately proportional 5,6 to sinQ for small 9, and may 
therefore be written 5 in the convenient form 

T d = -g(9)(h/2eSiSl)IS 2 x (Si x S 2 ), (1) 

where g(9) is a slowly varying function of 9. Note that the direction of 
a Gilbert or Landau-Lifshitz damping torque is also radial 9 . However, it is 
directed toward the precession axis of S2, while the current-induced drive torque 
Td is directed outwards provided I flows (Fig. 2a) from F\ toward F 2 , i.e., I > 0. 

For two electrons moving in different directions inside F 2 (Fig. 2b), the 
phases of s precession are equal at first, but become different at a given distance 
x from the N — F2 interface, as x increases from zero. Because of boundary 
conditions at the interface, the phase for an electron moving nearly parallel to 
the x axis changes with x at a slower rate than the phase for an electron going 
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in an oblique direction, assuming a spherical Fermi surface. Because of this 
dephasing effect, the torques exerted on S2 by different electrons tend to cancel 
5 for x larger than a certain characteristic length L = n/ (fcj — fcj.), where fcf, fcj. 
are the spin-up and spin-down wavenumbers in F2 (Fig. 2b). From the data of 
Ref. 4, and assuming a reasonable free-electron Fermi velocity of 1 x fO 6 m/s, 
we derive L ~ 9.5 nm and H sd ~ 1900 T for cobalt. The length L plays in 
the present problem the same role as the transverse spin relaxation time T 2 in 
nuclear magnetic resonance. 

Since regions of F2 with x > Lq do not contribute to the total torque on F2 
(except 8 for some periodic fluctuations), the current-induced torque per unit 
volume of F2 is inversely proportional 5 ' 6 ' 8 to the F2 thickness if if if > Lq. 
Electron scattering by solutes or phonons also reduces 6 the contribution of 
regions with large x to the torque, and has been observed to change s in the 
transmitted beam to a direction making a smaller angle with S 2 - 

III. SPIN CURRENT AND SPIN IMBALANCE 

In order to calculate the total drive torque fd on layer F 2 (Section II), one 
can introduce 5 the concept of angular-momentum current P(x), often called 
spin current. It is related to fd by fd — P(x = 0) — P(x = L%). I n practice, the 
second term is of no importance if Lf ^> L 0l f° r the same reasons (Section II) 
that regions with x > Lq do not contribute to the total torque. This approach 
has the advantage that the details of the exchange interaction between s and S2 
do not appear explicitly in the calculations. Also, it applies even to very large 
9 values. 

As an alternative approach, one can introduce 6 the difference Afi = ^ — /x j 
between the Fermi levels /ij, ji[ of spin-up and spin-down electrons. In Fig. 2c, 
we show the occupation numbers /j , f± versus energy e for spin- up and spin- 
down Fermi distributions at T ~ 0, if A/i < 0. Quantum theory 9 shows that 
the energy of a spin wave in F 2 consists of an integer number n m of energy 
quanta hu> (magnons), where oj is the S2 precession frequency. Note that n m 
is directly related to the average 9 angle, assumed to be small. The simplest 
collision process (Fig 3a), caused by s-d exchange and conserving angular mo- 
mentum, is 10 such that an electron flips its spin s from down to up while one 
magnon of energy huj and wavevector q is created. We show a complicated in- 
teraction vertex at the center of Fig. 3a, as a reminder that the simple Born 
approximation (Section I) is not valid for this process. Energy conservation 
gives ej — e-f = Tioj, where e-j-, are the final and initial electron energies in this 
so-called "spin-flip process" , also represented by an oblique arrow on Fig. 2c. 

Because of the Pauli exclusion principle, such processes must start from a 
filled state and end in an empty state. Fig. 2c shows that the only suitable 
initial states for these processes are those located in a region of energy thickness 
|A/i| —%uj = —(A/i + hw) near the top of the spin-down Fermi sea. Thus, the 
number of possible initial states and, in turn, the resulting drive torque are 
proportional to Afi + Uoj: 

f d = a sW (A M + M^#^. (2) 
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The term A/i is proportional to the current I, and thus is consistent with 
Eq.(l). But the new Hoj term represents 6 a surface damping torque indepen- 
dent of I, with a s (9) as the corresponding dimensionless Gilbert parameter. It 
is found that a s is inversely proportional to Lf. This surface damping is a 
consequence of the Pauli exclusion principle, and is not obtained if the spins 52 
are assumed classical and static as in the Slonczewski theory 5 . Although the 
present argument leading to Eq. (2) is based on a Fermi distribution with sharp 
Fermi level at T ~ and with Au + huj < (Fig. 2c), Eq. (2) turns out to be 
valid 6 even at high temperature fc^T » hu),Au where the Fermi level is very 
blurred, and even when A/i + Tilo > 0. 

The current I flowing from F\ to F% (Fig. 2a) may act 5,6 through a trans- 
lation of the Fermi surfaces in momentum space. In ferromagncts, the spin-up 
and spin-down currents, as well as the corresponding Fermi-surface displace- 
ments (Fig. 3b), are different, leading to non-zero P and A/i at a given point 
of the Fermi surface, with \Au\ < 1CP 4 eV . We saw in Section II that the im- 
portant electrons are those originating in F\ and moving toward F 2 . They are 
located on the k x > half of the Fermi surface, and indicated by a filling of dots 
on Fig. 3b. By summing over all such k x > states, one can calculate the total 
spin current T.P(x = 0) and the average A/i. Actually, only one component 
of P, here called P x2 , is important. It is found 11 that spin current and A/i 
remain in a constant ratio as the current I is varied, or as the ratio of spin-up 
to spin-down conductivities is varied: 

ZP x2 (x = 0) 

— — = const.. 3) 

An ■ sinB 

This shows that the Slonczewski theory 5 , based on T,P x2 , and the Berger 
theory 6 , based on A/i, are largely equivalent. The most significant difference 
between them is that, in their original form, the first one uses phenomcnological 
bulk damping only, and the second surface damping only. This difference has 
observable consequences, which will be described in Section V and Fig. 4b. 

There exists another kind of spin imbalance, associated with isotropic expan- 
sion or contraction of the Fermi surfaces (Fig. 3c), which goes back to Aronov 
and to Johnson and Silsbce 12 . It is the same at all points of the Fermi surface. 
We denote it by A/i, and it can reach 10~ 3 eV. Like A/z, it is related linearly 
to I, but the relation is more nonlocal. It exists only within a spin-diffusion 
length from interfaces. As before, only electrons on the k x > half of the Fermi 
surface are really active in creating torques, and they are shown by a filling 
of dots in Fig. 3c. Again, the total spin-current component and average spin 
imbalance, now denoted by TiP x2 and Au, can be calculated 11 . Again, we find 
SP Z 2/ A/i = const., and the constant has nearly the same value as before. 

Numerical estimates 11 , based on the assumption of one-dimensional cur- 
rent flow, suggest that T,P x2 and Au for Fermi-surface expansion/contraction 
are each 10 or 100 times larger than the corresponding quantities for transla- 
tion. Thus, the expansion/contraction mechanism is dominant if the flow is 
one-dimensional, i.e., if current leads attached to the multilayer sample do not 
flare out to a larger diameter too close to the multilayer itself. 

IV. CRITICAL CURRENT VERSUS FIELD 
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In Section II, we saw that the current-induced drive torque fd on layer F 2 
(Eq. 1) is directed outwards, i.e., tends to increase 0, if I > 0. On the other 
hand, bulk Gilbert damping, characterized by the Gilbert parameter a b , gives 
an inward torque: 

s 2 ds 2 £ 2 x(£ix£ 2 ) 

r b = -a b — X — = a fc w — . (4) 

o 2 o\o 2 

Strictly speaking, the second equality only holds for circular precession in 
a field H normal to layers. The physical origin of this bulk damping will be 
discussed briefly in Section VII. 

At small I or |A/z|, damping dominates, so that 9 actually decreases as time 
passes (Fig. 3d). However, if I is positive and exceeds a certain critical value I c , 
the drive torque overcomes damping, and 9 increases (Fig. 3e). Depending on 
the function g{9) in Eq. (1), 9 may stabilize at a constant value, or may go all the 
way to 180° and stay there. Detailed numerical as well as analytical calculations 
of dc current-induced precession and switching have been done 13 for various 
thin-film geometries and anisotropy cases. And even when the precessing spin 
distribution in F 2 becomes non- uniform 14 . 

The value of I c is obtained by writing t<j + f& = 0, where we take Td from Eq. 
(2) rather than Eq. (1) to include surface damping. Using also Eq. (4), and 
writing a s — C s /Lf an d A^i = —KI, where C s and K are positive constants, 
we obtain 

/ c = ^(l + ^)^. (5) 

In the case of H normal to layers, with circular S 2 precession, we have 9 a 
frequency lo = j(H z — M s ) for the uniform mode. Here, M s is the saturation 
magnetization in SI units. There is an additional 47r factor in front of M s if cgs 
units are used. Eq. (5) shows that I c itself is proportional to lo: 

/ c cx W = 7 (ff z -M s ). (6) 

Here, 7 is the gyromagnetic ratio. Tsoi et al. 15,18 measured I c on a multi- 
layer consisting of many alternating Co and Cu layers, normal to the field. A 
sharp silver point (Fig. 4a) contacting the sample was used to achieve large local 
current densities. Consistent with Eq. (6), they found I c to increase linearly 
with increasing H Zl though not quite with the predicted vertical intercept. 

The case 9 of in-plane field H with u = ^(H Z (H Z + M S )) 1 ^ 2 is more common, 
as it does not require large fields to saturate the sample, but more complicated. 
Since S 2 precession is now elliptical (Fig. 3f), the bulk damping torque ft, (Eq. 
4) and the two parts of fd vary through a precession cycle, in different manners. 
As a result, I c is not simply proportional to frequency w, but rather 16 linear in 
H z : 

I c kH z + (M s /2). (7) 

Again, an extra 4n factor should be put in front of M s if cgs units are used. 
The prediction of Eq. (7) is more or less consistent with I c measurements by 
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Katine et al. 16 and Albert et al. 17 on Co/Cu/Co multilayers in in-plane fields, 
where I c is found to increase linearly with increasing H z , with a positive vertical 
intercept. 

Tsoi et al. 18 have succeeded in measuring the frequency ui of current-induced 
spin oscillations. At fixed H z , they find that uj increases with increasing I above 
I c . This is not surprising: Since I c increases with increasing w (Eq. (6)), the 
only spin wave excited at / ~ I c is the one with lowest ui, i.e., the uniform 
mode with zero wavenumber q (Fig. 3g). At / > I c , other modes with larger 
uj and q may appear. Slonczewski 19 has even proposed that the excited modes 
at / ~ I c already have q > 0, due to the radiation of spin waves away from the 
point contact; this is rather consistent with the measured u) at / = I c in Tsoi's 
Fig. 3. 

To measure u>, Tsoi et al. adapted a technique well known to radio engineers 
under the name of "synchronous rectification" or of "direct conversion". The 
sample is exposed to a microwave of known and adjustable frequency oj ex (Fig. 
4a). While varying the dc current I through the point contact, they measured 
the dc ohmic voltage V between point and sample. Due to nonlinear effects, a 
glitch appears 18 in dV/dl whenever u ex is equal to the frequency u> of current- 
induced spin precession. 

V. CRITICAL CURRENT VERSUS LAYER THICKNESS 

According to Eq.(5), I c is a linear function of Fi thickness Lf , with positive 
intercept on the I c axis. This prediction is shown as the upper solid line in Fig. 
4b. Actually, it should hold for any direction of H. 

I c has been measured versus Lf by Albert et al. 20 , for in-plane field. They 
found I c to vary linearly with L 2 , but with no appreciable vertical intercept. If 
we were to ignore the surface-damping term Twj in Eq.(2), the constant first term 
would disappear from Eq. (5), bringing that equation into better agreement with 
the experimental findings above. This is shown as the lower solid line in Fig. 
4b. 

On the other hand, Urban et al. 21 have recently measured the Gilbert 
damping constant versus precessing-layer thickness in very thin Fe/Au/Fe/Au 
multilayers with single-crystal layers, using a transmission-resonance technique. 
They find a considerable contribution inversely proportional to thickness, in 
other words, a surface damping. Conceivably, the discrepancy between the two 
experiments may come from a non-magnetic "dead" atomic layer 4 reducing the 
effective Lf . 

VI. NEXI TORQUES IN VERY THIN LAYERS. 

We saw in Section II that the total current-induced torque on the spins S 2 of 
layer F 2 is usually in the radial direction t</ of the S 2 precession orbit of Fig. 2a. 
However, for very small F 2 thickness L 2 < L — 10 nm 7 the spin s of an electron 
does not have enough time to precess appreciably around S 2 while crossing F 2 
at the Fermi velocity (Fig 2b). Rather, s remains largely in its initial direction 
parallel or antiparallcl to S\. Therefore, the total torque is in the tangential 
direction f n (Fig. 2a), like the torque of a current-induced exchange field of 
fixed direction parallel or antiparallel to Si. This effect was first proposed by 
Hcide et al. 22 , who called it non-equilibrium exchange interaction (NEXI), and 
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expected it to exist at all layer thicknesses. 

We have calculated the magnitude of r„ as a function of if, by the same 
methods 8 as for t<j- The normalized torque per unit volume of F 2 is plotted 
versus kjqL^ in Fig. 5, for k^jk^ — 1.5 and fc^/fcjv = 1-0. Here, fc;v is the Fermi 
wavenumber in layer N. The corresponding plot for Td was in Fig. 2 of Ref. 8. 
For Lf < L , T n is found to be comparable to fd- However, as expected, f n 
vanishes at larger L% values, except for periodic oscillations around zero. The 
period of these oscillations is approximately 2L = 2n/ (fcj — fcjj on the L% scale. 

VII. ANISOTROPIC S-D EXCHANGE IN THE BULK 

We saw in Section I that ordinary s-d exchange is not active in the bulk. 
Anisotropic s-d exchange is 23 a combination of spin-orbit interaction and ex- 
change, and does not conserve total angular momentum. It is probably respon- 
sible for most of the Gilbert damping (Eq. (4)) observed in bulk and single- film 
magnetic metals. This theoretical expectation is confirmed by the observation 
24 of an increase of damping in Ni and Co at low temperature, of the kind 
predicted 23 on the basis of anisotropic s-d exchange. 

VIII. ANALOGY WITH LASER DIODE 

Because of a rather complete analogy with the semiconductor laser, we called 
6 SWASER (Spin Wave Amplification by Stimulated Emission of Radiation) the 
current-induced spin instability described in the present paper. Like the rate 
of creation of photons in a laser, the rate of creation dn m /dt for magnons is 6 
proportional to the number n m of magnons itself, and this is called stimulated 
emission. By opposition, spontaneous emission is independent of n m , and is 
probably negligible 25 for magnons whenever 9 ^> 0.1 °. The energy gap 7-ff z in 
the spin- wave spectrum plays in our theory a role like that of the semiconductor 
band gap in a laser. Also, the spin imbalance A/x is the analog of the difference 
u v ~ fi c between the Fermi levels 26 of the valence and conduction bands. The 
spin- wave damping proportional to Tilo (see Section III) also exists for light waves 
in a laser 26 . As a result, our condition 6 Au + Tilo = for the critical current 
has 26 a correspondent u v — [i c + Tilo = for the laser threshold. Finally, the 
N — F2 interface acts as a source of momentum 8 to bridge the momentum gap 
fc-f — fc| between spin-up and spin-down Fermi surfaces. Similarly, solute atoms 
and lattice pores are used 27 to bridge the momentum gap between valence and 
conduction bands in lasers made of an indirect-gap semiconductor such as Si. 
One difference, however, is that spin waves are far more nonlinear than light 
waves. 

The SWASER can be used 17 to reliably switch small magnetic elements, 
leading to possible applications for data storage. Also, it constitutes a very 
compact, tunable, microwave oscillator. 
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FIG. 1. a) Spin configuration assumed in the Born approximation, with s 
nearly static and antiparallel to H. b) More correct picture with s precessing 
rapidly around S 2 , and following it closely, c) Observation of s precession by 
an angle e around a fixed S, in Fe, Co, Ni films (Ref. 4). 

FIG. 2. a) Multilayer with magnetic layers Fi,F 2 and nonmagnetic N,N 2 . 
The oblique arrow represents an electron crossing from F\ to F 2 . b) Two elec- 
trons propagating inside F 2 in different directions have different s precession 
phases around S 2 at a given x > 0. c) Fermi distributions /|(e), /j.(e) for spin- 
up and spin-down electrons of energy e, with different Fermi levels at 
T ~ 0. The oblique arrow represents a spin- flip process where a magnon of 
energy Uco is emitted. 

FIG. 3. a) Spin-flip process where a magnon of energy hco and wavevector 
q is created, b) Current-induced translations of spin-up and spin-down Fermi 
surfaces, c) Current-induced expansion or contraction of spin-up and spin-down 
Fermi surfaces, d) Precession orbit of spin S 2 if I < I c . e) Orbit if I > I c . f) 
Orbit if I > I c and H is in-plane, g) Dispersion relation u(q) for spin waves in 
a field H z . Horizontal lines show the maximum oj value of current-induced spin 
waves. 

FIG. 4. a) Multilayer with many magnetic layers and a point contact, used by 
Tsoi et al. (Refs. 15 and 18). b) The upper solid line represents the predicted 
dependence of I c on F 2 thickness L 2 if both surface and bulk dampings are 
introduced. The lower solid line applies if only bulk damping is introduced. 
The crosses represent schematically the experimental data of Albert et al. (Ref. 
20). 

FIG. 5. Normalized tangential torque component r n per unit volume of F 2 , 
plotted versus normalized layer thickness L 2 , for fixed I, 9 and u. The point 
corresponding to L 2 = L is marked on the horizontal scale. 
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